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Abstract 

We show that supercritically charged black holes with NUT provide a new setting for traversable 
wormholes. This does not require exotic matter, a price being the Misner string singularities. Without 
assuming time periodicity to make Misner strings unobservable, we show that, contrary to expecta¬ 
tions, geodesics do not stop there. Moreover, since there is no central singularity the space-time turns 
out to be geodesically complete. Another unpleasant feature of spacetimes with NUTs is the presence 
of regions where the azimuthal angle ip becomes timelike, signalling the appearance of closed timelike 
curves (CTCs). We show that among them there are no closed timelike or null geodesics, so the freely 
falling observers should not encounter causality violations. Considering worldlines of charged particles, 
we find that, although these can become closed in the vicinity of the wormhole throat for large enough 
charge-to-mass ratio, the non-causal orbits are still disconnected from the distant zones. Integrating 
the geodesic equations completely, we demonstrate the existence of timelike and null geodesics con¬ 
necting two asymptotic regions of the wormhole, such that the tidal forces in the throat are reasonably 
small. We also discuss bounds on the NUT charge which follow from the Schwinger pair creation and 
ionization thresholds. 
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I. INTRODUCTION 


In this paper we would like to explore a novel aspect of non-vacuum space-times endowed 
with a Newman-Unti-Tamburino (NUT) parameter. Recall that the vacuum Taub solution was 
hrst presented as an anisotropic cosmological solution [1], and later rediscovered by Newman, 
Tamburino and Unti [2] as a static black hole (for more details see [3]). This latter solution 
has a unique event horizon of the Schwarzschild type, the corresponding internal metric being 
just the original Taub cosmology. Soon after Brill presented a corresponding pair for the 
Einstein-Maxwell system [4], its black hole face is commonly called the Reissner-Nordstrom- 
NUT (RN-NUT) solution. Depending on the relative values of parameters, it can have two 
horizons, one degenerate horizon, or no horizon. The cosmological Brill solution corresponds 
to the region between two horizons in the hrst case. Here we investigate the third case of 
no horizon, which surprisingly was not discussed in the literature so far. It corresponds to a 
Lorentzian wormhole connecting two asymptotically locally hat spaces. Because it is free of 
a central singularity, there is no place for a mass or charge, so that the Brill wormhole is a 
realization of Wheeler’s “charge without charge” and “mass without mass” [5]. 

An astonishing feature of this “NUT wormhole” is that it does not demand an exotic mat¬ 
ter violating the null energy condition (NEC), or the averaged NEC (ANEC), crucial for the 
existence of spherical [6-8] or non spherical [9] asymptotically hat wormholes. Indeed, it is 
supported by a Maxwell held which is not exotic. But the price for that is the presence of a 
metric singularity on the polar axis known as the Misner string, which is analogous to the Dirac 
string of the magnetic monopole in electrodynamics. To make the string unobservable Misner 
suggested to impose periodicity of the time coordinate [10], which entails quantization of the 
energy of matter helds [11, 12] similar to Dirac charge quantization. This, however, leads to 
violation of causality throughout the space-time, so it can hardly be considered as a physically 
acceptable condition. In the black hole case, another problem is that analytic continuation 
cannot be consistently carried out through both horizons, so that the resulting spacetime is 
geodesically incomplete. In the wormhole case there are no horizons and analytic continuation 
just reduces to the extension of the semiaxis or the radial variable r to the whole axis. So in the 
case of the Misner interpretation the wormhole RN-NUT spacetime is geodesically complete. 

In an alternative interpretation of the Misner string suggested by Bonnor [13] (see also [14]), 
the time periodicity condition is abandoned and the polar axis is treated as a kind of topological 
defect with its proper matter source. We will assume this viewpoint here. Since the polar axis 
is now an unremovable singularity of spacetime, it was expected so far that geodesics will stop 
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here. Recently we have shown [15] that it is not true for the pure Taub-NUT metric. Here we 
extend this proof to the Brill charged solution as well. Our proof is independent on whether 
one assumes the Misner or Bonnor interpretation of the space-time. In the wormhole case the 
space-time structure is particularly simple, and since there is no singularity other than the 
polar axis this means that the Brill spacetime is geodesically complete. 

Also, in the Bonnor interpretation there is a region around the Misner string where closed 
timelike curves (CTCs) exist. Recently we have shown that in the case of the Taub-NUT 
spacetime these CTCs are not geodesic for a suitable choice of a parameter hxing the position of 
the Misner string. Thus a freely falling observer will not be confronted with causality violation. 
Here we extend the proof to the Brill solution too. Since we deal with charged wormholes, it is 
natural to consider also the worldlines of charged particles. We show that a circular worldline 
lying at the throat of a massless magnetically charged wormhole can be causality violating for 
large enough charge-to-mass ratio and a certain sign of the charge. At the same time, we show 
that all the wordlines extending to large distances from the throat are causal. Together with 
the proof of geodesic completeness, this makes us believe that NUT wormholes are free from 
traditional objections against solutions with NUTs. 

Passing to the analysis of the physical features of the solutions, we hrst demonstrate the 
existence of timelike and null geodesics connecting two asymptotic zones. Then we investigate 
the tidal forces in the vicinity of the throat and show that these may be reasonably small for 
a large NUT charge. More restrictive bounds on this parameter are obtained by demanding 
that the electric field in the throat is lower than the Schwinger pair creation or the matter 
destruction characteristic helds. We conclude with some suggestions for further theoretical 
work. 

II. THE BRILL SOLUTION 

The solution of the Einstein-Maxwell system of equations found by Brill in 1964 [4] soon 
after the discovery of Newman, Tamburino and Unti [2] reads 

ds'^ = —f{dt — 2n(cos 6 + C) dcpY + f~^dr‘^ + (r^ -|- n^) [dO"^ + sin^ Od^p ^), 

A = ^(dt — 2n(cosd + C) d(p), (2.1) 

with 

(r — m)^ + gr-h p(r^ — n^)/2n 

J= + ’ = r2 + n2 ’ 

/l2 2 I 2 2 2\ 

[b = q + p — m — n ). 
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The metric depends only on the combination of the electric (g) and magnetic {p) 

charges. For e = 0 (6^ = —(m^ + n^)) the solution reduces to the Taub-NUT solution. For n = 0 
(6^ = — m^), it reduces to the dyonic Reissner-Nordstrom solution. It is particularly simple 

in the massless case m = 0 for the charges satisfying = 2n^ (6^ = n^), then the gravistatic 
potential /(r) = 1, and the solution looks like a “pure” NUT. Finally, for + p^ = + ri^ 

(6 = 0), the solution is a special case of the Israel-Wilson-Perjes (IWP) solution: 

ds^ = —f [dt + A.dS^ + f~^dx ‘^, R = <F (^dt + A.dSj , 
f~^ = 1+ 2acosa + , V A ^ = 2 sin aVa, 

<h = —-[siu/S-f 2(Tsin(Q;-|-/9)-|-sin(2Q;-|-/9)], (2.3) 

^ SlKl 

where (j{x) a harmonic function. 

The value of the parameter C introduced above can be modihed by a “large” coordinate 
transformation 

t^t + 2n{C-C')p. (2.4) 


Such a transformation will generically lead to an intrinsically different spacetime. Thus, (2.1) 
actually dehnes a one-parameter family of spacetimes, to which we will refer collectively as “the 
Brill spacetime”. We will keep this parameter C free until section 6 . 

In the case of a non-vanishing NUT charge, n 7 ^ 0, the Brill solution (2.1) is not singular at 
r = 0 , as can be checked by computing the quadratic curvature invariants, 

4e^ 






R^^upaR 


{A + ’ 

8 e^ 


(2.5) 


+ 


48 


{(m^ — n^) [r® — 15n^r"^ -|- 15nV^ — n® 


(A + n^ ' {A + Ay 

—2mr [(e^ — Gn^A — (lOe^ — 20n^)n^r^ -|- (5e^ — 6A)A] 
-fe^[(e^ — 10n^)r^ — 2(3e^ — 10n^)n^r^ -f (e^ — 2A)A]] . 


( 2 . 6 ) 


The Kretschmann scalar (2.6) (where the first and second terms are respectively the Ricci 
square and Weyl square contributions), reduces to that of the Reissner-Nordstrom solution [16] 
for A = 0, and to that of the Schwarzschild-NUT solution [17] for = 0. In view of its non¬ 
singularity, the Brill solution can be considered as a regularization of the Reissner-Nordstrom 
solution. For 6 ^ < 0 it has, just as the RN solution, two horizons. For A = 0, it has, just as 
the extreme RN solution, a double horizon. However for 6 ^ > 0, contrary to the RN solution, 
it is not singular, but has the (Lorentzian) wormhole topology, the coordinate r varying in the 
whole real axis, with two asymptotic regions r = ±cxd. As r > 0 decreases, 2 -spheres r = 
constant shrink until a minimal sphere of area 47rn^ (the wormhole neck) for r = 0, and then 
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expand as r < 0 continues to decrease. Because of the absence of a point singularity, there 
is no source for the various (gravi-)electric and (gravi-)magnetic fluxes, and thus this solution 
realizes the Wheeler program [5] of mass without mass, charge without charge, etc. (for a more 
recent discussion see [18]. Note also that (as in the Taub-NUT case e = 0), the mass parameter 
m is not positive dehnite, as the reflection r ^ —r changes its sign and that of the electric 
charge q. 

The price to pay for this regularization of the RN point singularity is the introduction of 
the Misner string singularity. The “Misner string” (which coincides with the Dirac string for 
p 7 ^ 0) consists of two disconnected inhnite components, the North string piercing the North 
poles 0 = 0 of all the spheres r = constant (r G R) and its counterpart the South string for 
6 * = TT, except in the special gauges C = =f 1, where only one component (South or North) 
is singular. The Misner string singularity can be transformed away altogether if the time 
coordinate t is periodically identihed with period Smr. However there is also a price to pay 
for this. First, there are closed timelike curves (CTC) everywhere in all the Brill spacetimes. 
Second, in the case 6 ^ < 0, the Kruskal continuation cannot be consistently carried out with this 
periodical identihcation, so that the black-hole Brill spacetimes are not geodesically complete 
[3, 19], i.e. the Misner string singularity has been traded for a singular horizon. 

So let us keep the usual real time axis, thereby retaining also the Misner string. In this case 
it has been argued [19] (in the case of the Taub-NUT black hole spacetime, but exactly the 
same argument could be formulated in the general Brill case) that geodesics must terminate 
at the Misner string singularity, just because this is a metric singularity. Actually, geodesic 
motion does not depend on the metric itself, but only on the connections, so a metric singularity 
does not necessarily imply geodesic termination. Geodesics do indeed terminate at a conical 
singularity (the analogy made in [19]), but not at the metric singularity r = 0 of the plane with 
metric df' = dr"^ + r'^dO'^l The investigation of geodesic motion in the metric (2.1), carried out 
in the next section, will show 1) that all geodesics which hit the Misner string cross it smoothy, 
so that all the Brill spacetimes with 7 ^ 0 are geodesically complete; 2) that if furthermore 

> 0 the resulting Lorentzian wormholes are traversable. 
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III. GEODESICS: ANGULAR AND TEMPORAL MOTION 


The partially integrated geodesic equations for the metric (2.1) are [20] 

/(r ){i — 2n(cos 6 + C) (p) = E , 

. LM 

^ (^2 _|_ ^ 2 ^ gj ]^2 Q ’ 

[(r^ + v?)9] ■ = (r^ + n^) sin 9 cos 9 pp' — 2nE sin 6^ 0 , 
f{r)~^{P — E‘^) + (r^ + rP){9‘^ + siiP 9 (p‘^) = e , 

where ' = d/dr, with r an affine parameter, 

Lz{9) = Jz — 2nE cos 9 , (3.5) 

and E and Jz are the constants of the motion associated with the cyclic variables t and ip (see 
Appendix A). In the following we assume without loss of generality that for timelike and null 
geodesics i? > 0, thereby dehning the orientation of time, and n > 0, the sign of Jz remaining 
arbitrary. Eq. (3.4) is ds^ = edr"^, with e = —1, 0 or 1 for timelike, null or spacelike geodesics, 
respectively. 

The analysis of geodesic motion in the Brill metric parallels that of a charged particle in 
the held of a magnetic monopole [20]. Similarly to the case of the Taub-NUT metric, the 
metric (2.1) admits four Killing vectors, one {Kq) generating time translations, and three {Kp 
generating the rotation group 0(3). Due to this spherical symmetry, the angular components 
of the geodesic equations can be hrst integrated by 

L + S = J, (3.6) 

where J is a constant vector, the total angular momentum, which is the sum of an orbital 
angular momentum L and a “spin” angular momentum S. The components of L and S are 
given in [20], where it is shown that 

L = (r^ + n^) f A f , S = 2nEf (3.7) 

where 

r = (sin 6^ cos 9 ?, sin6'sin(/), cos6') (3.8) 

is a unit vector normal to the two-sphere. 

It follows from the orthogonality of L and S that 

J-f = 2nE. (3.9) 


(3.1) 

(3.2) 

(3.3) 

(3.4) 
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In the magnetic monopole case, such a first integral means that the the trajectory of the charged 
particle lies on the surface of a cone with axis J originating from the magnetic monopole source 
r = 0. However the Taub-NUT or Brill gravitational held has no source, i.e. no apex for 
the “cone”. Actually, the content of (3.9) is that the geodesic intersects all the two-spheres of 
radius r on the same small circle, or parallel, C with polar axis J. Squaring (3.6) leads to 

= P + 4:n^E'^ , (3.10) 

which can be rewritten as 

+ Pf[e‘^+ siPep] = i\ (3.11) 

with P = P — 4pE‘^ [P = P). Insertion into (3.4) leads to the effective radial equation 

r‘^ + f{r) J = E^, (3.12) 

which is identical to the equation for radial motion in the equatorial plane for the metric ( 2 . 1 ) 
without the term —2n cos 6dip, i.e. for 6 ^ > 0 motion in a static spherically symmetric wormhole 
geometry of the Ellis type [21] (the m = 0 symmetrical Ellis wormhole [22] for /(r) = 1). 

A. Angular motion 

Knowing (in principle) the solution to (3.12), one can insert it in (3.11), where 0 can 
be eliminated by using (3.2), to yield a hrst-order differential equation for 0{t). Then the 
solution can be used to obtain ip^r) and f(r) by integrating (3.2) and (3.1). The Misner string 
singularities 6 * = 0 , vr do not seem to play any role in the complete integration thus carried out. 
To check this, we follow [19] and replace the affine parameter r by the new variable A dehned 
by 

dr = [P + P) d\ , (3.13) 

which increases monotonously with r. Putting 

^ = cos 9 , (3-14) 

the differential equation for 6{X) reads 

+ + Jf). (3.15) 

Assuming 7 ^ 0 (J^ = 0 implies from (3.10) E = 0 and I = 0), Eq. (3.15) is solved (up to 
an additive constant to A) by [19] 


cos 9 = J ^ [2nEJz -f IJ± cos(JA)] , 


(3.16) 



where 

J| = J2 _ j2 _ 

Because the constants of the motion occurring in (3.16) are related by the two Pythagorean 
decompositions of (3.10) and (3.17), it will be useful to trade them for the two angles rj and 
tjj dehned by 


2nE = J cos rj, / = J sin 77 , Jz = J cos ip , J± = J sin tjj , 


(3.18) 


with 0 < ?7 < 7 r /2 and <ip <tt. Then (3.16) reads 


cos 6 * = COS'?/'cos ?7 + sin'^sin? 7 cos(JA). 


Eq. (3.15) has two turning points 6± such that 


cos 6± = J {2nEJz ± IJ±) = cos('^ =F rj) 


(3.19) 


(3.20) 


These can be rewritten as 


9^ = 


Ip — 7] \i Ip >7] 
—Ip + 7] if Ip <7j 


9 . = 


pJ + 7] if Ip <71 — 7] 

—Ip — 7] + 2% if Ip >7T — 7] 


(3.21) 


It follows that the trajectory crosses periodically the Misner string, cos 9± = ±1 only if 

Ip = 7] {Jz = 2nE) or ip = n — rj {J^ = —2nE). (3.22) 

The only geodesics which can cross both components of the Misner string are those with t] = 
7 r /2 {E = Jz = 0 ), leading to i = <p = 0] according to (3.12), in the stationary sector 
(/(^) > 0) these can only be spacelike geodesics. The trajectory can also stay on the Misner 
string component 6 ^ = 0 or vr if (3.22)) is satished with 77 = 0 {2nE = J). 

The differential equation (3.2) for cp can be rewritten as 


dp 


Jz — 2nE 


+ 


Jz + 2nE 


d\ 2 [l —cos 6 *(A) l + cos 6 '(A)J ’ 
with cos 6 *(A) given by (3.16). This is solved by [19] 


(3.23) 


p — Po = arctan 


arctan 


cosip — COST] J\ 

---- tan — 

1 — cos('^ — 77) 2 

cos'^ +cos 77 JA 

I- 7~i -T tan — 

[1 + cos('d’ — 77) 2 


with pq an integration constant. This can be simplihed to 

2 sin 77 tan(JA/ 2 ) 


p — Pq = arctan 


[sin(77 — Ip) — sin(77 + ip) tan^{JX/2) 


(3.24) 


(3.25) 


9 



For trajectories crossing the North Misner string, with 'ip = r] = 


ip — (pi = arctan 


^cos ?7 tan 





2nE), this reduces to 


(3.26) 


with Pi = pq — sgn(tan( JA/ 2 )) 7 r/ 2 . A similar formula applies in the case of the South Misner 
string, with 7] replaced by tt — 77 and J\ replaced by JA — tt (note that according to (3.16) 
the North Misner string is crossed for A = 2kTi/ J, while the South Misner string is crossed for 
A = {2k + l) 7 r/J, k integer). In the case e.g. of the North Misner string, this gives on account 
of (3.16), 


J. 


cos{p - Pi) = — tan - 


(3.27) 


consistent with (3.9) (the choice 9 ?! = 0 in (3.27) corresponds to the choice J = (Jx,0, J^) in 
(3.9)). Clearly the Misner string is completely transparent to the geodesic motion! 

When the parameter A varies over a period, e.g. A G [—n/J, tt/ J], the argument of the hrst 
or second arctan in (3.24) varies from —cx) to +cx) for =|= 2nE > 0, and from +cx) to —cx) for 
Jz =F 2nE < 0. It is identically zero for Jz =F 2nE = 0. Accordingly, the variation of p over a 
period is 


A 79 = TT [sgn( — 2nE) + sgn( + 2nE)] . 


(3.28) 


This means that for J/ > Ain?E‘^ (I^T’I = 27r) the parallel C circles the North-South polar axis, 
i.e. the Misner string. For J/ < AmpE"^, (1^7^I = 0) C does not circle the Misner string. And 
for Jz = ±2nE (| A (^| = tt), C goes through the North or South pole, as discussed above. 

The two turning points (3.21) coincide if either ip = 0, or r] = 0: 

aP) Tp = 0 or TT. The constant vector J is aligned with the 2 : axis, on which the parallel C is 
centered: 


6 = rj or TT — rj, p — pii = ±J\. (3.29) 

As discussed in [20], the solution (3.16), (3.25) can always be rotated to this case, the Misner 
string being then rotated away from the 2 axis. 

b) ?7 = 0. The orbital angular momentum vanishes, / = 0, and the motion is purely radial. 


6 = 7p , p- pii = {). 


(3.30) 


This will be further discussed in section 4. 


B. Temporal motion 

After transforming to the geodesic variable r by (3.13), the solution to (3.1) can be written 
as the sum 

t{\) = tr{\) + t0{\) , (3.31) 


10 



where the radial and angular contributions to t(A) solve the equations 


dtr ^r‘^ + n^ 

dte 2n{cos 6 + C){Jz —‘InE cos 6) 

dX sin^ 9 


(3.32) 

(3.33) 


The solution to equation (3.32) depends on the solution of the equation for radial motion (3.12). 
We consider here equation (3.33). This can be rewritten as 

dte 


dX 


= An E + n 


{C + 1)(J, - 2nE) _^{C- 1)(J, + 2nE) 


1 — cos 9{X) 


1 + cos6'(A) 


with cos6'(A) given by (3.16). The explicit solution to equation (3.34) is [19] 
te{X) = te{X) = An^EX + 2n{C + 1) arctan 
+ 2n{C — 1) arctan 


cos Ip — COST! JX 

. ---- tan — 

[1 — COS('^ — 7 ]) 2 


cos'^ + cos ?7 JX 

-—-- tan — 

1 + cos(y’ — Tj) 2 


(3.34) 


(3.35) 


in the interval —tt/J < A < tt/J. The resulting variation of te over a period 27i/J of A is 

AnE 


Ate = 27in 


J 


+ {C + l)sgn( - 2nE) + {C - l)sgn( + 2nE) 


(3.36) 


so that the solution for generic values of A is 


te{X) — te{X) + 



Sttu^E 

J 



(3.37) 


For parallels C which do not circle the Misner string (J^ < An‘^E‘^), Ate is negative and inde¬ 
pendent of the value of C, Ate = —47rn(l — cost]). In the limit when such parallels contract to 
a point, the motion becomes purely radial {2nE = J) and the angular contribution te vanishes. 

Let us recall here that, while the solution for the spatial components of the geodesic (the 
orbit) can be rotated to a frame in which the axis points along the total angular momentum 
vector J, this is not the case for the time coordinate t{X). The correct equation for the time 
evolution in this rotated frame is given (for the choice C* = 0) in [20], Eq. (24a). 


IV. GEODESICS: THE COMPLETE ORBITS 

A. The radial potential 

The radial equation (3.12) can be written in the familiar form 

(^) +U{r) = E^ + e, (4.1) 
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FIG. 1: Radial potential corresponding to n = 
1, m = .4, / = 2.8, b = .3 with two maxima and 
two minima. The rightest one at r = 16.23395, 
U = —0.2383e —1, not seen in this scale, is shown 
separately in Fig. 2. 



FIG. 2: The Newtonian potential well region of 
the potential Fig. 1 with higher resolution. 


with the effective potential 


f/(r) = £ 


2mr + /^[(r — m)^ + 6 ^] 






(4.2) 


with Its derivative generically can presented as U' = Pi{r )/where 

Pi{r) is a polynomial of fourth order in r, which, depending on the parameter values, may have 
four real roots or two real roots. In the cases m = 0 or e = 0, P 4 (r) degenerates to a polynomial 
of third order with three or one real roots. 

In the case of timelike geodesics {e = —1), for m 7 ^ 0 the potential is asymmetric under 
reflection r —r. We will assume m > 0 which corresponds to a positive mass from the point 
of view of an observer at r = +00 and negative mass for an observer at r = — 00 . In the generic 
case one maximum is at negative r while at r > 0 there are two minima with a local maximum 
between, or only one minimum. A typical potential with four extrema is shown in Figs. 1-2. 
The case with only two extrema is illustrated in Fig. 3. Remarkably, the potential for zero 
orbital momentum / = 0, generically has two extrema: a minimum at r = r_|_ and a maximum 
at r = r_ (see Fig. 4), 

— 2n? ± \Je^ — 

2m 

corresponding to stable or unstable equilibrium positions. The potential for the massless worm- 
hole m = 0 is symmetric under reflection r —r, and so has an extremum at r = 0, and two 


(4.3) 
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FIG. 3: The potential with two extrema for non¬ 
zero orbital momentum (n = m = / = 1, 6 = .3). 



FIG. 5: The generic potential for zero mass m = 
0 with three extrema and n = l,/ = l,6=.2 



FIG. 4: The potential for zero orbital momentum 
I = 0 and n = 1, m = .4,6 = .5. For I = 0 the 
extrema can be found analytically by solving a 
quadratic equation. 



FIG. 6: Massless wormhole potential with one 
maximum (m = 0,1 = 0.1, n = 1, 6 = 1.3). 


other extrema at r = ±ro: 

2 {2h‘^ — + in? ih'^ — v?) 

^ r?-h^- P ’ 

provided 0 < r? — I? <? and ? < n?{n? — I?)/{21? — r?) (if 46^ > r?/2). 

The potential for nnll geodesics e = 0 is positive dehnite (Figs. 8-10). It has either one, 
or three extrema. In the massless case m = 0 there is an extremnm at r = 0, and two other 
extrema at r = ±ri with 



r‘l = r? — 21? , (4.5) 

provided r? > 26^. In this case there is a stable circnlar photon orbit at the wormhole throat 
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FIG. 7: Massless wormhole potential with one 
minimum (m = 0, / = .1, n = 1, 6 = 0.1. 



FIG. 9: Massless wormhole potential for null 
geodesics with one maximum (m = 0,1 = 
0.1, n = 1, 6 = 1.3). 


/ ^ 

/ 1.5 

/ 0.5 - 

-8 -6 -4 -2 0 2 4 6 8 

r 

FIG. 8: Generic potential for null geodesics with 
three extrema (n = 1, m = 1/2, I = 2, b = .2). 



FIG. 10: Massless wormhole potential for null 
geodesics with three extrema (m = 0, I = .1, n = 
l,b = 0.1). 



r = 0, with areal radius n, and unstable photon orbits at r = ±ri. The purely radial motion 
{I = 0) of the photons is free. 

By following the procedure of [19], exact solutions for r(A) can be obtained in terms of 
Weirstrass elliptic functions. Here we will only discuss the qualitative radial motion for timelike 
or null geodesics sampling the wormhole (6^ > 0) geometry, and consider special limiting cases 
of interest. We shall work in the rotated coordinate frame in which the z axis is along J, so 
that the polar angle 6 is constant. 
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B. Timelike geodesics (e = — 1) 


1. Scattering states 

The effective potential is for m > 0 negative for large positive r and positive for large negative 
r, and so presents at least a maximum and a minimum, with > ''"max, and possibly also 
a secondary maximum and a secondary minimum. A test particle coming from inhnity with 
> 1 will first be accelerated in the potential well around rmm, before either being reflected 
(scattered) by the potential barrier or, if its energy is high enough, being transmitted through 
the wormhole neck r = r^ax fo the second asymptotic sheet r ^ —oo. 

Combining this radial motion with the angular motion previously discussed leads to a tra¬ 
jectory which spirals around the polar axis J at a constant angular velocity J with respect to 
A. The dehnition (3.13) of A together with (4.1) leads to the relation 


d\ = 


(r2 + n2)[E2-l-C(r)]V2 ' 


Putting — 1, it follows that the geodesic range of A is given by the hnite integrals 

f°° dr 

(r^ + n^W-U{r)YI^ 

in the case of reflexion (with tq the turning point, 17(ro) = /3^), or 

dr 

AAtrans = J (^2 + ^2) [^2 _ t;(^)] 1/2 

in the case of transmission (if = L Anax )- Dehning as usual the impact parameter by 

a = 1/Id, the cross-section for transmission is 

2 

O'trans = = — - . (4.9) 

^max 

In the limiting intermediate case (d'^ = (d^, the integral (4.7) with the lower bound rmax is 
logarithmically divergent, and the trajectory, attracted by the unstable circular orbit r = r^ax, 
spirals indehnitely towards the wormhole neck. In the present paper, we will only evaluate 
these integrals for two special cases of interest: the non-relativistic limit, and the small-NUT 


The non-relativistic limit is dehned by m = 0 and 6^ = (so that /(r) = 1) and 1. 

In this case, 

7/(r) =/^/(r^ 4-n^) (4-10) 


leads to r‘l = a? — vd, with 


A 1 2 7°° dr 

d Xo [(r2 + n2)(r2-r2)]i/2 ~ T \a 


(4.11) 
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where K is the complete elliptic integral of the hrst kind. In the frame '^ = 0, the corresponding 
total angular variation is 




K ( cot 9 


(4.12) 


where we have used i? ~ 1 for small f3. If / ~ 2ntan6' is held hxed in the non-relativistic limit 
(corresponding to large impact parameters a), (4.12) reduces (using i^(0) = 7r/2) to the hnite 
result Ai^refi — tt/ sin0. If on the other hand the impact parameter a is held hxed (of the order 
of n), then the cone angle 9 goes to zero, and the test particle spirals around the polar axis 
many times before being scattered back to inhnity. 

Similarly, in the case of transmission through the wormhole {a < ao = n), 


A A + r-o pc. 


= -^K 


P Jo [(r^ + n2)(r2 + — a2)]i/2 /jn 


(4.13) 


Again, if 1 the cone angle 0 ~ /3a/2n is very small, and the test particle makes a hnite 
but large number of turns 

(4.14) 

27r n/3 \n) 

along its path through the wormhole. 

Remarkably, the number of turns remains large, 

^ i ( 4 . 15 ) 

in the limit of a purely radial trajectory, / = 0 (a = 0). This fact, which results from the 
conservation of total angular momentum (3.6) (J = S' in the case of purely radial motion, 
L = 0), is of course not possible for a classical point test particle, but makes sense if one 
considers rather a small extended solid test body. For instance, (4.15) means that a tennis ball, 
of radius small before the wormhole throat radius n, will spin around, with an angular velocity 


d^p 2nE 
dr 


(4.16) 


which is maximum at the wormhole throat r = 0, making turns while going straight through 
the wormhole from r = +cx)tor = —cx). 

Let us emphasize that this uniform rotation of an extended test body in free radial fall does 
not depend on the direction of this radial motion, which generically does not coincide with 
that of the Misner string. From (3.19) and (3.25), the angular motion along a parallel C is 
a uniform rotation, with proper angular velocity J/(r^ + n^). In the limit of a purely radial 
motion, I —)■ 0, the parallel C contracts to a point 9 = constant, ip =constant, while the uniform 
rotation velocity retains the finite value 2nE/{r‘^ + n^). 
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The small-NUT limit is motivated by the non-observation of NUT charge in nature. It is 
dehned by v? tv? and 6^, so that ~ e^. In this limit the effective potential (4.2) 

can be approximated by 


U{r) ~ 


-2mr^ + (e^ -|- — 2mPr + 

(^2 _|_ 


and its derivative by 


U'{r) ~ 


2r[mr^ — (e^ -|- P)r‘^ -|- Sml'^r — 2e^P 

(^2 _|_ ^ 2^3 


(4.17) 


(4.18) 


Obviously r = 0 is a maximum of U{r), and an absolute maximum because U(0) — U{r) > 0 
in the small NUT limit: 


/32 = ~ U(0) ~ 


2/2 


eH 




7^ 0) 1 


^ (/ = o). 


(4.19) 




We hrst consider the case of reflexion (/9^ < Pq), which we will discuss only in the small-velocity 
limit (/9 —)■ 0). In this limit I = jSa can be neglected (purely radial motion), so that 

2mr — 




u ^-U ^ 


r2 _|_ ^2 


(4.20) 


with vq ~ e^/2m. The evaluation of the integral (4.7) yields a small, but non-vanishing 


27r 


2n 

e 


(4.21) 


Transmission through the wormhole neck occurs if /3^ > /3q, i.e. in the extreme-relativistic 
case. The cone angle is given in this case by tan0 ~ a/2n < n/2e. In the extreme-relativistic 
limit we can approximate 




(j,2 _|_ ^2y 


(3^ 





(^2 _|_ ^ 



(4.22) 


The evaluation of the integral (4.8) yields 


^V^trans — 


4n 




-.K 


2ea 


ea 


71^ + ea 


(4.23) 


For very small impact parameters (corresponding to hnite anglar momenta I = /Sa) a n^/e, 
(4.23) reduces to the hnite result Ai/^trans — 27r. Again, this remains valid in the limit of a purely 
radial trajectory, so that a small extended test body falling straight through the wormhole will 
at the same time spin around (at a maximum proper angular velocity dip/dr ~ 2/9/n), making 
exactly one turn while going from r = -|-cx) to r = —oo. 
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2. Bound states 


Bound states in the wormhole Brill geometry have some unusual features with respect to the 
typical black hole situation. The generic potential (Fig. 1) for massive particles has two poten¬ 
tial wells: one is at negative U in the far region with small binding energies (non-relativistic 
Newtonian orbits) shown in Fig. 2. Another, relativistic potential well closer to the throat, 
extends to positive energies, being separated from escape by a potential barrier (Fig. 1) with 
larger binding energies. One has therefore oscillating bound orbits between the two turning 
points of two types: Newtonian, and relativistic. At the corresponding minima one has two 
stable circular orbits. If two roots of the equation U\r) = 0 are complex, the potential has a 
simpler shape with only two extrema. Fig. 3, the minimum (potential well) being in the r > 0 
region, and the maximum at some r < 0 . As it is seen from this hgure, the potential well is 
deep and relativistic, the Newtonian well being absent in this case. 

For zero orbital momentum / = 0, the potential reduces to the simple form 


!/(.) = '’A 1 


(4.24) 


with a minimum and a maximum given by Eq. (4.3). One can observe that the the shape of 
the potential for I = 0 (Fig. 3) is qualitatively similar to that for / ^ 0 with two extrema (Fig. 
4). Indeed for / = 0 we still have a non-zero total angular momentum J = 2nE. The bound 
orbits are oscillating between two turning points with relativistic velocities. The test particle 


can remain at rest at the minimum at r = r_|_ provided its squared energy 



is positive, i.e. if > 0. However for n 7 ^ 0, a small test body at the equilibrium position 
r = tq will spin around with angular velocity Uq = 2nEo/(rQ -|- n^). In the small-NUT limit, 
these / = 0 equilibrium values reduce to 


^ b 2nbm^ 

ro ~ , Eo ~ , Wo --^ • 

me e° 


(4.26) 


The potential is particularly simple and reflection-symmetric in the case of the massless 
wormhole m = 0. Depending on the parameter values, there could be three different cases: a 
potential well centered at the throat extending to positive energy values and separated from 
the scattering region by potential barriers (Fig. 5), a pure potential barrier centered at the 
throat (Fig. 6 ) and a purely negative potential well at the throat (Fig. 7). In the first and 
the last cases one thus has stable bound orbits oscillating around the throat, including circular 
ones located exactly at the throat. 
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Another simple case is the near-extreme limit m?. The effective potential (4.2) has its 
absolute minimum at r = m (the extreme RN horizon radius) for 6^ = 0 (e^ = m? -|- n^). So 
for 6^ small and positive, there are stable circular orbits at r = m -|- 0(6^), with energy Ei and 
proper angular velocity ui = dip/dr 


E 


2 

I 




(4.27) 


C. Null geodesics (e = 0) 


The effective potential for null geodesics 


U{r) 


/^[(r — mY + 6^] 
(^2 _|_ ’ 


(4.28) 


is shown in Figs. 8-10. The discussion of the scattering of light parallels that for timelike 
geodesics. Of special interest is the scattering angle A^^refi = ■^AArefl, where A Area is given by 
(4.7) with P = E, and = {4rP -|- aP)E‘^. This can be computed analytically in the special 
case m = 0, with the result 


Av^refl 


Av^refl 



4 

T) 

for > -rp, (4.29) 

tP 

for rY > 2Y and A(rP — Y) < < —r , 

¥ 


with 

a^ = aYa^ + A{Y-nY]. rg = . (4.30) 

In the non-relativistic case 6^ = rP, this reduces for a > n to 

= 2^1 + A K (-) = P-K (i cot o) , (4.31) 

V ¥ V a / sm 6* \ 2 J 

where tan^ >1/2, while light rays with small impact parameters a < n (tan6' < 1/2) end up 
in the second asymptotic sheet (r —oo) with 

Av^trans = (2 tan 6^) . (4.32) 

COSO* 


Let us also note the existence, in the near-extreme limit rrP, of circular null orbits at 
r czi m (as in the case of massive particles) with energy 

P 


Ef ~ 6 — . 

t „A 


(4.33) 
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V. CHARGED PARTICLE MOTION 


Consider now the non-geodesic motion of a charged particle with mass rric and charge Qc- It 
is convenient to work with the covariant Lorentz equations of motion: 

, (5.1) 

where k, = qc/mc- Since has only g, = a = {t, ip) components depending only on r, 9, one 
has for yU = a 

Fa.x’^ = -Aa , (5.2) 

while for the remaining g = i = [r, 6) components 

F^^x'^ = d^AaxF (5.3) 

The equations of motion thus split into two groups. The first group of equations do not contain 
the metric derivative term and combine to a total derivative 

{x'"gau + nAa)'= Q , (5.4) 

which can be integrated as before by introducing two integrals of motion E, 

f{r){i — 2n{cos9 + C) p) — K^{r) = E , (5.5) 

x'^g^iy + kA^ = Jz + ‘^nC , (5.6) 

the second one being equivalent to (3.2). In the second group it is enough to consider only the 
9 component, taking the normalization condition (3.4) as a constraint. The i = 9 component 
reads 

1 . 1 

[(r^ 71^)9] ■ = -gabpx'^x^ + = gtippiip + + nA^^e p. (5.7) 

Eliminating t from (5.5), we hnd that the K-terms cancel, and we are left precisely with the 
equation 

[(r^ -I- rA)9 ]■ = (r^ -I- rA) sin 9 cos 9 p^ — 2nE sm9 p , (5.8) 

identical to (3.3). So the angular equations of motion are the same as in the case of a neutral 
particle, leading to orbits which he on a “cone” of half-angle arctan(2ni7//) (with I the orbital 
angular momentum), while in the equation for the time evolution (3.1) and in the radial equation 
(3.12) the constant E must be replaced by E -|- K^{r) leading to 

+ /W + ll = (E + K^r))^ , (5.9) 

- 1 - 
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where we took e = —1 for massive, charged particles. This can be rewritten as 

r P 1 

r^ + hh(r) = 0, where W{r) = f{r) — - 5 -+ 1 — (E + K^{r)P . (5.10) 

\_r^ + J 

The quantity W (r) however is not convenient for the role of the effective radial potential since it 
depends quadratically on E and contains unphysical negative energies. To correctly introduce 
the radial potential we have to present this equation in the factorized form 

= {E-V+{r)){E, (5.11) 

where 

r±(r) = --.t ± + (5.12) 

SO that for a particle at rest r = 0 the correct branch E = V+. Moreover, the quantity 
E + K^{r) is equal to the kinetic energy, it has therefore to be strictly positive (recall that 
/ > 0 everywhere for the wormhole): 

E + K^{r) = i r'^ + f{r) —^ + 1 >0- (5.13) 

W _j_ j^Z 

It follows that E > V- always, so the particle’s motion corresponds to > l + . Also, the 
conditions for circular orbits W = 0 = W are implied by the more physical conditions = 
Vl = 0 which do not contain unphysical negative energies. So we can consider V+ as a correct 
radial potential. 

To compare with the results of [23] for the motion of charged particles in the dyonic Reissner- 
Nordstrom spacetime, one must hrst gauge transform our potential <l>(r) to a potential vanishing 
at inhnity ‘ho(^")- This transformation and the associated shift in the integration constant E 
are 

$(r) = $o(r) + |-, E = Eo-'^. (5.14) 

Taking the NUT charge to zero, we hnd that the orbit of a charged particle in the dyonic 
Reissner-Nordstrom spacetime lies on a cone of half-angle arctan(/T;p//), in agreement with the 
results of [23]. 

^ Recall that in special relativity the total energy of a charge in the potential y (note that $ is dehned with an 
opposite sign to y) satisfying {E — ey)^ = m'^ + p‘^ has to be solved with the sign plus: E = eip+ -f p^. 
In our case the mass is absorbed by the affine parameter on the worldline, so to = 1. 
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FIG. 11: Radial potential for charged particle 
motion in the wormhole with m = 0, (? = 0, = 

2n^ in the case of positive k for different angular 
momenta l‘^ {k = \/2, = 19, 12, 6, 0 decreas¬ 

ing from top to bottom). 


FIG. 12: Radial potential for negative k (= —y/2) 
and the same values of P from top to bottom. 


A. Motion in the case / = 1 


We restrict to the particularly simple case of the massless wormhole with only magnetic and 
NUT charges related so tha /(r) = 1, 


m = 0, q = 0 , p = \/2n , 


leading to 


14 = 


2 2 
K r — n 


+ t 1 + 




(5.15) 


(5.16) 


i/2 ' V ' 

The potential V+ is symmetrical under reflection r —>■ —r. For k > 0 it is maximum at the 
wormhole throat r = 0 and monotonically decreases with growing |r| (Fig. 11). The maximal 
value depends on the orbital momentum as follows: 




(5.17) 


For K < 0 the shape of the potential is more diverse, as shown in Fig. 12. In the region I: 


— <2V2\k\, 


(5.18) 


the potential has a local minimum at r = 0. With growing P the local minimum goes up and 
in the region II; 


2V2\k\ < — < 44 + 24142(1 + 24), 


(5.19) 
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local maxima develop on both sides of the throat, with the positions 


r± = ±n 


SK^in?{in? + P) — P 

P - 


(5.20) 


Finally, in the region III, 

]2 , . 

^ + 2\K\^/2(l + 2K^)j , (5.21) 

the two local maxima at r = r± disappear, while the local minimum at r = 0 turns into an 
absolute maximum. 

With the potentials shown in Fig. 11, 12 one deduces the following characterization of the 
orbits. For k > 0 one has orbits traversing the wormhole throat for E > Vq corresponding to 
the top of the potential barrier in Fig. 11. The orbits with E < Vq correspond to scattering on 
the wormhole. For k < 0 in the parameter region III the situation is the same. In region II one 
has in addition bound states between the turning points inside the potential well between the 
local maxima. The same energies correspond to scattering orbits reflected on the exterior sides 
of the barriers. The region I contains bound orbits inside the well and traversing orbits for the 
energies above the wells. These can be analyzed in analogy with Sec. 4, we leave it for future 
work. 

The circular orbits at the throat r = 0 (with areal radius n) exist for both signs of n. These 
will be interesting for analysis of non-causality in the next section. They are unstable for 
positive K and negative k in the region III, but stable for negative k in the regions I and II. The 
corresponding gauged energy E = h+(0) -|- k/\/ 2 is positive for positive k, but it is negative for 
negative k if 


^<2k^ 


corresponding to the kinetic energy 

E + K^{r±) = 

and positive total energy 

^ 1^1 ( 
a/2 


{\k\ > 1/V2). 

(5.22) 

circular orbits exist 

in region II at r = r± 

P 

(5.23) 

2^/2\K\n^ ’ 

2r?2\ 


(5.24) 


VI. CAUSALITY 


The ADM form of the metric(2.1) is 


ds^ = 


/(r^ -I- rP) sin^ 9 


dP + f ^dP + {P + n^)d9^ + [ dcp + 


2nf{cos9 + C) 


dt ) , (6.1) 
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with 


S(r, 9) = (r^ + n^) sin^ 9 — 4n^/(cos 9 + C)'^ . 


( 6 . 2 ) 


For /(r) >0 (which is always the case for the wormhole Brill solution with 6^ > 0), S becomes 
negative, and closed timelike curves (CTCs) appear, in a neighborhood of the Misner string 
given by S(r, 0) < 0. The surface S(r, 0) = 0 bounding this CTC neighborhood is a causal 
singularity of the spacetime, where the signature of the spacetime changes from (—1-++) outside 
to (+ + H—) inside. This singularity is, just as the Misner string itself, completely transparent 
to geodesic motion. Nevertheless, the occurrence of CTCs in a spacetime is usually considered 
to violate causality [3, 24]. An observer traveling around snch a CTC wonld eventnally retnrn to 
his original spacetime position after a hnite proper time lapse, thus opening the possibility for 
time travel. However, unless this observer is freely falling, snch a CTC travel would necessarily 
involve accelerations generated e.g. by rocket engines. One can argue that the back-reaction 
of these matter accelerations on the spacetime geometry will deform it in snch a way that 
chronology will nltimately be preserved. If this reasoning is correct, in vacnnm gravity causality 
violation can only occur in spacetimes with closed timelike geodesics (CTCs), or possibly closed 
nnll geodesics (CNGs). However we deal here with a charged solntion, so we must also consider 
the worldlines of charged particles. 


A. No closed timelike or null geodesics 

We now show that there are no closed timelike or nnll geodesics in the Brill spacetimes with 
\C\ < 1. Closed geodesics (or self-intersecting geodesics) occur if after a hnite lapse of affine 
parameter r or of Mino time A, all the coordinates take again the same valnes (modnlo 27r 
for the azimnthal angle). So the perimeter of the geodesic mnst be an integer mnltiple of the 
A-period 2ti/ J of the angnlar motion. From (3.36), dnring this period, 

- l) (6.3) 

for \C\ < 1. Also, from (3.32) and (3.12), 

A = _ e(T^ + n=)] > E-'f - en ^], (6.4) 

d\ f{r) 

leading to 

Ac A ^ , (6.5) 

over the same period. Adding the two together, we obtain 

277 - 

At = Atr + Ate A ^ [J — 2.nE — en^/ J] . (6.6) 


Ate > 47rn 


/ 2nE 

w 
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For e = —1 this is clearly positive definite. For e = 0, this can vanish only for 2nE = J {I = 0). 
But in this case Ate > 0, while dtr/d\, and thus also Atr, is positive definite. Thus, for \C\ < 1 
all timelike or null geodesics are causal (future directed). 

The situation for |C| > 1 is less obvious. We shall only discuss the case of null geodesics. 
As we shall show, 1) there are always CNGs in any given Brill spacetime if \C\ is large enough; 
2) for any given value of \C\ > 1, there are Brill spacetimes with CNGs. 

First, we observe that for C* > 1, 


Ate = Ann 


2nE 

~r 



(6.7) 


for orbits which circle the South Misner string {J^ < —2nE). The same relation (6.7) holds 
for C* < —1 in the case of orbits which circle the North Misner string > 2nE). We wish to 
show that, if the orbit is closed, this negative angular contribution {2nE < J) can be exactly 
balanced by the positive radial contribution, leading to At = 0 and thus to a closed geodesic. 
Depending on the values of the model parameters, two situations can lead to closed orbits with 
e = 0. Either U{r) has a local minimum f/min > 0 between two local maxima, and there are 
bound states in the resulting potential well, leading to closed orbits if the periods of the radial 
and angular motion are commensurate (i.e. for a dense set of values of the energy). Or U{r) 
has only a maximum, corresponding to an unstable circular orbit. Let us discuss these two 
cases separately. 

In the first case, from (3.32) and (3.12), 


dtr ;2 ^ i2 ^ 


d\ U{r) 

It then follows that, over an angular period. 


< 




( 6 . 8 ) 


At < 47rn 


1 + 




2nE 


|C| 


< 47rn 


1 + 


f 


An^Urr 


|C| 


(6.9) 


J 

vanishing for 101 = 1 + /^/4n^f/min (which depends only on the wormhole parameters (m, n, b)). 
For higher values of |0|, the bound state null geodesics will be past directed. 

In the case of a circular orbit of radius tq and energy Eq = 77(ro) (circling the North or 
South Misner string as above). 


dtr 

dX Eq 


P - An^E^ 


En 


leading to 


2 n 

At = —[J-2nE,\C\], 


( 6 . 10 ) 


( 6 . 11 ) 
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which vanishes for 


1^1 = 


J 


2nEn 


= 1 + 






1/2 


4n2f/(ro )) 


( 6 . 12 ) 


(value depending only on the wormhole parameters). So again in this case, for any given set of 
Brill parameters {m,n,b), there will be CNGs for large enough \C\. 

We now formulate the converse, and stronger, proposition: that for any given \C\ > 1, one 
can find a set {m,n,b) such that there are CNGs. Presumably this also extends to CTGs. 
Actually this proposition is easy to prove in the small NUT limit. There is an unstable circular 
orbit at r ~ 0, of energy Eq ~ for I ^ 0. From the above, there is a family of 

corresponding closed null geodesics for 

,2 \ 1/2 


n 


which proves the proposition. 


|G| ~ 1 + rc , 


(6.13) 


B. Closed charged worldlines 


Since we deal with a charged solution, it is natural to also consider the worldlines of charged 
particles. In this case we have 


27r 

Ac > , 

^ J 


(6.14) 


which is positive as in the geodesic case by virtue of (5.13). The angular contribution Atg still 
satisfies Eq. (6.3) (restricting to \C\ < 1) and can be negative. Gonsider circular orbits at the 
wormhole throat r = 0 in the case /(r) = 1, g = 0, p = \/2n. From (5.13), Atr = 27rn7/J, 
with 7 = \//2 + n^, so that the total At is equal to 

27in 


At = 


J 


(7 + AnE — 2J) , 


(6.15) 


with J = A/ 4 n^E 2 + 7^ — n^. Denoting A = (7 + AnE^ — AJ"^ and using E = V+{fd) given by 
(5.17) we obtain 


A N/ A 2\/2K,n , ^ _ 

A = 5(7 - 7_)(7 - 7_+), 7± =--— |1±\/1 


Since 7 must satisfy 7 > 1, we find that A is non-positive for 

9 


5 

^2 


(6.16) 


K < 


■ 4 ^ 2 ’ 


7- < 7 < 7+ • 


(6.17) 


Note that the upper bound (6.17) is smaller than A\/2\K\n/b, leading to P/n'^ < 32 k^/25, which 
is much lower than the upper bound (5.19) of region II, so that the circular worldlines r = 0 
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satisfying (6.17) are stable. We expect that nearby bound worldlines in the potential well above 
r = 0 will similarly violate causality. However, such worldlines are conhned to a narrow region 
around the wormhole throat, so that these causality violations are actually unobservable, for 
an observer initially living at a large distance ri from the wormhole. 

Now let us argue that all the worldlines which can be followed by such an observer returning 
to his starting point are causal. The negative angular contribution to the total round trip time 
lapse Aftot is given by the contribution (6.3) during one period multiplied by a hnite number 
of periods proportional to 


AA 


refl — 



dr 

(r2 + n^)[{E - H + (r))(E - H_(r))]V2- 


( 6 . 18 ) 


This will be easily balanced by the positive radial contribution which is proportional to the 
distance between the turning point tq and ri, so that the worldline will remain causal in the 
large. The observer travelling in a charged rocket and returning to his starting point will not 
meet his younger self. 

However this reasoning fails if, for a given value of the orbital angular momentum /, the 
traveller’s energy is just a bit below the maximum of V+(r). The number of A-periods (6.18) 
diverges logarithmically for tq = rmax so that, for tq close to rmax, the charged particle makes 
many turns around the wormhole throat before returning towards inhnity. Then, the angular 
time lapse could be very large (and negative), and the worldline could possibly self-intersect at 
multiple radii r = r^, if the attractor unstable worldline at r = r^ax was closed. But we have 
shown above that the wordlines at r = 0 can be closed only if they are stable. In the region H 
parameter range (5.19), there are unstable circular orbits at r = r±. During a period, the total 
At is 


At 


27r 

T 


(r^ + n^){E + fi;*h(r±)) -|- An^E — 2nVAn‘^E‘^ + P 


(6.19) 


where r± is given by (5.20), E + K<h(r±) by (5.23) and E by (5.24). To see whether the square 
root term may dominate we calculate the difference of the squared quantities 


{r^+n^){E+K^{r±))+An^ Ef —An^ {An^ E^+r) = 


(H - 


+ 2{y + An‘^){l^ — Sk^u'^) , 


( 6 . 20 ) 

which is positive by virtue of Eq. (5.19), so that these unstable circular wordlines are causal. 
So, while this does not constitute a rigorous proof, we believe that there are no closed charged 
wordlines extending to spacelike inhnity. The case of wormhole parameters more generic than 
those (5.15) considered here should similarly be investigated. 
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VII. PHYSICAL ASPECTS 


A. Tidal accelerations and traversability 


In order to be traversable, macroscopic Lorentzian wormholes should not only be geodesically 
complete, with a class of timelike or null geodesics connecting the two asymptotic regions, but 
they should also be such that the tidal gravitational forces exerted on an observer travelling 
through remain reasonably small [6]. We will discuss here only the case of radial geodesics 
(/ = 0). In the non-relativistic case (/(r) = 1), the radial velocity is constant, so the longitudinal 
tidal force will vanish. However we have noted previously that an extended object falling 
radially through the wormhole will rotate at a uniform proper angular velocity 

2nE 


u{r) = 


(7.1) 


r2 _|_ ^2 

which should generate a centrifugal acceleration a(r). In flat space, this would be given in 
cylindrical coordinates (p, z) by 

7 .. ^ 4n^ 


aflat = p = uj [z) p 2 Z 


where we have assumed for simplicity that the direction of motion is along the axis, approx¬ 
imated r zz z (the radius of the object should be small before the wormhole throat radius n), 
and E zz 1 (non-relativistic limit). Taking as in [6] the reference values pmax = 2m for the 
radius of the object and Omax = 9 = 9.8 ms“^ (Earth gravity), we arrive at the rough estimate 
at the throat r = 0 

n A n _ 

(7.3) 


(7.2) 


-( 0 ) 4 

P 'Tl Pmax 


i^max r —2 

~ ~ 5 S 


leading to a NUT charge n ~ Is, corresponding to a wormhole throat radius n ~ 3.10®m (or in 
mass units, n ~ IO^Mq). 

To improve on this flat-space estimate, we must compute the tidal acceleration between 
two neighboring points of the freely falling test body with purely spatial separation 5x in the 
comoving frame [6], i.e. such that 

= 0 . (7.4) 

This tidal acceleration is given by the geodesic deviation equation [25] 


5x^ = -—6x^ = —R^p,y„u'^5x^u° 
dr^ 


(7.5) 


Before proceeding, we note that in the non-relativistic case the length scale is set by the sole 
NUT parameter n, so that the Riemann tensor components at the throat, and hence the relative 
tidal accelerations, must be of order as in our flat-spacet estimate. 



In the case of radial geodesics the four-velocity has only two non-vanishing compo¬ 

nents and u^, so that the condition (7.4) reduces to 


[goo^i + go’i^ip) + gu^r ^ 0 , 


(7,6) 


which may be used to eliminate St from the right-hand side of (7.5). Taking into account the 
non-vanishing values of the Riemann tensor components listed in Appendix A, we thus obtain 
for the longitudinal and transverse accelerations 


Sr = - [R\trStu\^ + R\rtSr{u^f + R\^rS^u\^] 
(/ — l)(3r^ — n^) -I- 4mr 

(r^ -|- n^)2 ’ 

Se = - [R^etiuy + R\er{u^?] SO 
r{r — m) — (r^ — n^)/ 

= (r2 + n2)2 


(7.7) 


(7.8) 


and a similar equation for S(p. In the non-relativistic case / = 1 (implying m = 0), the 
longitudinal acceleration (7.7) vanishes identically, while the transverse acceleration takes the 
value at the wormhole throat r = 0 


50(0) = 4 50. (7.9) 

n2 

This is precisely one fourth of the flat space estimate (7.3), meaning that our traversability 
estimate for the NUT charge should be halved, n ~ 0.5 s ~ 1.5 x 10® m. At present we have no 
simple explanation for this factor 1/4. 

For generic wormhole parameters, the longitudinal and transverse accelerations at the worm- 
hole throat take the values 

= -—^Sr, 

g2 — 

50(0) = - —se. (7.10) 

The transverse (centrifugal) acceleration is always positive > rn^), while the longitudinal 

acceleration can have either sign. In the small-NUT limit, the two accelerations become large 
(of order e^/n^) and opposite in sign, so that they work together to compress the test body 
along the direction of infall. 


B. Electromagnetic fields 

In the present case, the Lorentzian wormhole geometry is generated by electric and/or mag¬ 
netic helds which will become large near the wormhole throat. For simplicity we discuss only 
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the case of purely electric or purely magnetic monopole fields. In Gaussian units these helds 


are given by 


E = 
B = 


01 


= <1: 


r‘i _ .^2 


(r^ + 

1 r‘^ — -n? 

r^23 =P-, 


(P = 0) , 

(g = 0) , 


(7.11) 


with > V? + vY for NUT wormholes. We hrst consider a purely electric held, in 

which case q > n, so that the maximal held is Umax > A hrst constraint to be satished 
for the validity of the classical approximation which we have used is that the Schwinger pair 
creation should be unobservable. In hat space, the Schwinger critical held in h = c = 1 units 
reads Usch = he/^O; where pe, are the mass and the charge of the electron. The ratio can be 
put into the following form 

M© mpi f mpi'' ^ 


Usch n M© V /ie / ’ 

where mpi is the Planck mass. Inserting here ttipi/Mq = 1.1 x 10“^®, mpi/pe = 2.4 x 10^^, cq = 
13-17-1/2 obtains 

^ = 0.54x10®^. (7.13) 

A-Sch n 

Thus in the purely electric case the classical picture is valid if n ^ 5 x 10^ M© (which is of 
the order of the geometric radius of the Sun i?© ~ 7 x 10® m). This lower bound for n, which 
should be improved by considering pair creation in the wormhole spacetime metric is only six 
times larger than our traversability estimate. 

We should also consider the effect of large electric helds on infalling test bodies. While 
ordinary matter is electrically neutral, so that its motion is not affected by electromagnetic 
helds, it is made up of charged atomic nuclei and electrons, which are sensitive to these helds. 
For a spaceship to go through the wormhole without damage, the electric held at the throat 
should be smaller than the ionization threshold. One can take as the relevant electric held the 
Coulomb held of the electron at the Bohr radius 


E,, = 


60 

r'i ’ 
' B 


rp = 


l^eel 


(7.14) 


One hnds in the same units: 


= el = 137"® = 0.4 X 10"®. (7.15) 

Esch 

Therefore the threshold for traversability of an electric wormhole by an unmanned spaceship 
is n ^ 10^^ M© ~ 10^® m ~ 0.1 It-yr. Traversability by a manned spaceship would impose 
even more drastic constraints. For instance, molecular air is ionized by a constant electric held 
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E 3 X 10® V m“^ ~ to avoid this the throat radius n should be at least equal to 10® 

It-yr, of the order of the radius of our Galaxy, which seems totally unrealistic. 

In the case of a purely magnetic wormhole no pair creation occurs, but above the Schwinger 
threshold one can expect phenomena of superstrong magnetic held such as in certain neutron 
stars. For the wormhole radius value n ~ 1.5 x 10® m derived in the previous subsection, the 
wormhole-throat held values would be 5(0) ~ 3.3 x 10^® T. For comparison, human exposure to 
a constant magnetic held of up to 5 ~ 8 T can be tolerated [26]. Taking the latter (magnetic) 
value as reference would lead to the lower limit for the wormhole radius n > 6 x 10 ^^ m ~ 60 It- 
yr, corresponding to a huge NUT charge, larger than our previous tidal force estimate by 
nine orders of magnitude. However, traversability of a purely magnetic NUT wormhole by an 
unmanned spaceship would seem to be feasible for not too large wormhole radii. 


VIII. OUTLOOK 

In this paper we have suggested a novel type of wormholes withont exotic matter. They are 
supported by matter sources satisfying the NEC condition and evade the NEC-violation theorem 
due to the presence of a non-diagonal (t, 99 ) component of the stress tensor, associated with a 
NUT charge. The price to pay for this is that the spacetime is not strictly asymptotically hat, 
but only locally asymptotically hat, and presents a Misner string singularity. The particular 
realization described here is supported by a Maxwell held and corresponds to a supercritically 
charged RN-NUT solution without a central singularity or horizons. Extending the results 
of [15], we have analyzed in detail the geodesic motion and shown that the Misner string is 
transparent for geodesic motion without periodic identihcation of time, so that the spacetime 
is geodesically complete. This result is independent of the actual form of the gravitational 
potential /(r), and holds also for the black-hole and extreme black-hole RN-NUT solutions. A 
curious ehect of the NUT charge is that small freely falling test bodies are endowed with a proper 
angular momentum or spin, which is aligned with the direction of infall, and is independent of 
that of the Misner string. 

Even without periodic identihcation of time, the RN-NUT spacetimes contain regions with 
closed time-like curves. We have shown, however, that for some subfamily of these spacetimes 
there are no closed timelike or null geodesics, so that freely falling observers should not en¬ 
counter causality violations. On the other hand, the analysis of the motion of charged test 
particles has revealed that there are wordlines which can be closed or become past directed 
if their charge-to-mass ratio is large enough. This was found, however, for orbits threading 
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the throat and not reaching asymptotic regions of space-time, so the existence of observable 
manifestations of non-causality still remains open. 

We have also investigated under what conditions these RN-NUT wormholes are macroscop- 
ically traversable, and shown that the gravitational tidal forces at the throat could be kept 
reasonably small for a reasonably small NUT charge. However, the electromagnetic forces act¬ 
ing on ordinary matter under the same conditions would be extremely large, and could be kept 
under control only in the case of huge NUT charges. 

This work could be extended in several directions. First, our analysis was purely classical. 
One should explore propagation of waves and quantum implications of NUT wormholes, as well 
as their stability [27]. Also, similar NUT wormhole solutions could presumably be constructed 
in other gravitating held theories. One such solution was found recently in the Einstein-Skyrme 
theory [28] after our paper was completed. Five-dimensional vacuum gravity is known to admit 
NUTless Lorentzian wormhole solutions [29], which are however not traversable [30]. We suggest 
that the possible existence of geodesically complete and traversable NUTty wormhole solutions 
to this theory be investigated. Other hve-dimensional theories, such as hve-dimensional minimal 
supergravity (which can be consistently truncated to the four-dimensional Einstein-Maxwell 
theory [31]), should also presumably admit such solutions. 

Moreover, by superposing NUT-anti-NUT conhgurations one could try to construct asymp¬ 
totically hat wormhole spacetimes. In the extreme black-hole case (6^ = 0), one can construct 
stationary Israel-Wilson-Perjes NUT-anti-NUT superpositions. These are singular if time is 
periodically identihed [32], but can be regular (with only a Misner string connecting the two 
event horizons) if it is not. From the present analysis, they should be geodesically complete. 
It would be interesting to extend this to the construction of slowly orbiting NUT-anti-NUT 
wormhole dynamical solutions in the near-extreme case, i.e. for small b^> 0. These should be 
asymptotically hat, with only a non-contractible Misner string loop threading the two worm- 
hole throats. If the parameters could be arranged so that this two-wormhole conhguration 
with two regions at spacelike inhnity presented a high degree of symmetry, we speculate that it 
might be possible to identify these two asymptotically hat regions [33-35] to yield a traversable 
Wheeler-Misner wormhole [36] connecting two distant parts of a spacetime with only one region 
at inhnity. 

We conclude with some astrophysical speculations. Recently it was suggested that wormholes 
could exist in the galactic centers [37, 38] or in the outer regions of halos [39, 40]. The option 
to have a wormhole without exotic matter seems attractive: its mass and NUT charge could 
be reasonably of the scale of the mass concentration in the galactic centers (10® — 10 ® ATq), in 
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which case it would not be destroyed by quantum effects. However, crucial for this wormhole is 
the supercritical electric or magnetic charge, whose origin in astrophysical conditions is unclear 
for the moment. Nevertheless we think that further study of new realizations of gravitational 
fields with NUTs may be useful for future astrophysical applications. 
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Appendix A: Symmetries 


The metric (2.1) has the following four Killing vectors A"(a), cl = x,y,z,t: 


K(x) 

^(y) 

Ki.) 

Kit) 


2n(l + C cos 6) cos ip 
sin 9 

2n(l + C cos 6) sin p 
sin 6* 


■dt — sin pde — cos p cot , 
dt + cos pde — sin p cot 6d^ , 


+ 2nCdt , 




t ? 


(A.l) 


the first three forming the algebra so{3). These formulas extend previously known ones [10-12] 
to arbitrary G-parameter. Note that the presence of (9t-terms in the so{3) subalgebra reflects 
the necessity of a compensating time-shift while performing spatial rotations. The algebra 
so(3) can not be integrated to the group S'0(3), but leads to the unitary representation of the 
group SU{2), provided periodic identification of time according to Misner is performed. This 
is similar to the case of the magnetic monopole [41], for further references see [11]. 

The vector potential one-form A = ^{dt — 2n{cos6 + C) dp) is also symmetric under the 
action of these four isometries, the Lie derivatives are zero for all a: 


Ck, = 0 = Ki 


(a) 


(a) 




These symmetries generate four conserved quantities 


= K'^a)iKg^,v + kA^) , 


(A.2) 
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which read explicitly: 

/(a;) = 2nE{sin 6)~^ cos If — sin 6— cot 6 cos (f srn^ 6(p + 2nE cos 6^ , 

/(y) = 2nE{sin9)~^ sin f + cosOO"^ p^ — coiO sin f (^p^ sin^ 9if + 2nE cos9^ , (A.3) 

I{z) = p^ sin^ 9if + 2nE cos9 , (A.4) 

lit) = E. 

They do not depend on the vector potential. 


Appendix B: Riemann tensor 


Here we present the Newman-Penrose invariant projections of the Ricci and Weyl tensor for 
the metric (2.1). Choosing the null tetrad as 


^ iw 

"" “ 71 (77’ “'A’ ° 



(A.l) 


= 


we have one non-zero projection of the Ricci tensor 




(A.2) 


and one projection of the Weyl tensor: 

4/2 = C^uXTl^nffh^rE = 


{rn? + b‘^){r‘^ — v?) + n‘^r{3m — 2r) 

(^2 _|_ 

. 2nr{'w?‘ -|- b‘^){r‘^ + n^) + nr'^{r — 3m) + n^{r — m) 

+ *- T ^.—577 -• (A-3) 

These do not depend on the polar angles. The components of the Riemann tensor can be read 
off easily. We give explicitly the ones used in Sec. 7: 

,(/ — !) (3r^ — n^) -|- 4mr 




ttr 


= -/- 


(^2 _|_ ^2^2 


R\r^ = R\rt = -2n (cos 9 + C)f 


(/ — 1) (3r^ — rR) + Amr 


(^2 _|_ ^2y 


R 


tte 


= R^tt^ = -f 


rjd _ r)ip _ 

it rrB rnp 


r {r — m) — {R — rR) f 
{R + rR)^ 
r {r — m) — {R — rR) f 
f{R + n‘^)‘^ 


(A.4) 
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